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' Within the framework of periodic asymmetric Anderson model for Kondo isoulators an 

£N) ■ effective singlet-triplet Hamiltonian with indirect antiferromagnetic f — f exchange interaction 

is introduced which allows to study analytically the dynamic magnetic susceptibilities X/(k, 
of /-electrons. The approach allows to describe the three-level spin excitation spectrum with 
a specific dispersion in YbBi 2 . Distinctive feature of the consideration is the introduction of 
5^ , small radius singlet and triplet collective f — d excitations which at movement on a lattice form 

low - and high-energy spin bands. 
PACS: 71.28. +d,71.27.+a, 75.30.Mb 
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a. 

t3 : 1 Introduction 

a . 

Kondo-insulators (KI) represent a special class of strongly correlated systems. They demonstrate 
such phenomenon as intermediate valency, a narrow gap (or a pseudo- gap) in a spectrum of carriers 
of an order of 10 meV which opens at low temperatures, an unusual spectrum of spin fluctuations 
(with a spin gap), a singlet ground state and a set of other interesting properties (reviews [T], [2], [3], 
[I], [S], |6J, [7]). The theory of KI is based on the periodic Anderson model (PAM) which describes 
the intrasite hybridization V between d— electrons and localised / -electrons with intrasite Coulomb 
/ — / repulsion Uf. The intersite interaction answers to the hopping Hamiltonian t for d— electrons. 
Considering low-temperature spin fluctuations, we will discuss such a characteristic KI as YbB 12 . At 
T < 40i\~, close to the spin gap edge the compound YbB u demonstrates three dispersive excitations 
Ml, M2, h M3 with average energies 15, 20 h 40 meV (see [8J and references in [8]). The relaxation 
of low-energy excitation Ml is suppressed, that allows to consider it as resonance excitation. Peak 
Ml has a dispersion with a minimum, and peak M3 with a maximum at L point which is the 
03 ' aniferromagnetic (AFM) point for the face-centered cubic YbBi 2 lattice. 

Let's give the known accepted approaches to treatment of a magnetic susceptibility x(k, a;) of 
YbB\ 2 . The particle-hole symmetric limit of the PAM with /— level filling nj = 1 is often used 
[lj,[9| (however this mode is unjust in the mixed valence regime (MVR) |TU]) . Sometimes strong 
/ — / correlation is considered in the slave boson mean-field approximation and x(k, u>) calculations 
within the framework of PAM use random phase approximation with introduction of crystalline 
electric field split for 4/— states [TT], [TJ], [IB]. And at last, it was also proposed an alternative 
model which is explicitly based on dielectric singlet ground state of YbB 12 for which it is taken 
that Yb -site f 14 , f^d 1 configurations are nearly degenerate [H]. Although all the treatments are 
devoted to the microscopic origin and fine structure of the spin gap in YbB\ 2) in the majority of 
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cases x(k, uj) considerations have a phenomenological and qualitative character, are generally based 
on the single-impurity Anderson model, so that AFM correlations between localized electrons arise 
in higher order perturbation-theory approximations in hybridization and hopping interactions. In 
particular, as far as we know, the explicit form for x(k, to) -dispersion of the lower spin excitation 
branches was studied only in [13], [15]. However there was no explanation of high-energy peak M3. 

In present paper it will be investigated an asymmetrical PAM variant (APAM) with two electrons 
per site in MVR rif > 1, assming that it is realised in YbB 12 - 

In case of MVR the magnetic system actually does not have small parameter. The close situation 
is realised in doped CuOi planes in HTSC cuprates. where the indirect d—d intersite AFM exchange 
is of the same order as p-holes interaction with the spin subsystem. In this case it is known, that for 
the adequate account of strong correlations for the low-energy excitations first place is necessary to 
solve a problem for small radius cluster, next to construct low-energy Zhang- Rice polaron (collective 
excitation of small radius) [IE] (see also [17]) and only in the end to consider such a polaron movement 
on a lattice. 

Similarly, for the accepted model it will be shown, that observable three-level spin excitation 
character can be described if from the very beginning the problem is considered on the basis of one- 
and two- site f — d collective triplet-singlet excitations (small radius excitations). Consideration 
differs from the previous themes, that x(k, uj) is defined exactly by movement of these excitations 
against a singlet ground state (movement is determined by the indirect AFM /-/ exchange). 

Below we will be interested only in a susceptibilityx/(k, uj) of /- electrons assuming that at low 
temperatures it is main for full x(k, a;). The effective singlet-triplet Hamiltonian H s _ t is itroduced 
for two site-cluster analysis and X/(k, uj) calculations. It is substantiated that H s _ t describes the 
basic properties of APAM. 

Initial background of our approach are close to the phenomenological x(k, uj) consideration in 
[T4] . however as a whole the approach is more realistic and it allows not only to explain Ml, M2, 
M3 peaks, but also to find their spectrum. 

2 Effective Hamiltonian and two-site cluster. 

As shown in [15] for the description of KI spin subsystem with rif > 1 it is convenient to use the 
effective Hamiltonian Hj which turns out from PAM Hamiltonian by a replacement of a hopping 
turm t on AFM / — / indirect exchange interaction J on the nearest neighbors. 
In conventional notation Hj Hamiltonian has the form 

Hj = H + J; J = - J ^2 S/,rS/ ir+g , H = ^2 Ho,n', (1) 

r,g n 
H , n = [^2£ffn,afn,a + Uff+ + fn,+fn,-fn,-] + ^ 6 dd+^dg^ + V n ] 

CT (7 

V n = Vj2(fn,X,« + dtJn,«); = ±, 

a 

g is a nearest neighbor vector. Here we neglect orbital degeneracy. V and t - amplitudes of V u t 
interactions; £/ and Ed are energies of /—and d— levels, below we put Ed = 0; J - a constant of AFM 
/ — / exchange. MVR with rif > 1 is characterized by the following energy parametres relations: 
Uf + £ f ~ V > and V < U f ~ -e f . 

Let us discuss a full set of two-electron eigen states of one-site Hamiltonian H r (pQ) in order to 
validate an introduced below effective singlet-triplet Hamiltonian H s _ t (it will differ from jl]) by of 
an intrasite term). 



2 



At T < — Ef and the accepted parametres it is possible at once to omit a high -energy state state 
2rf r = d^ + d^_\0) (|0)- site vacuum) with two particles on d - level. 

There are one singlet r and three triplet states with one electron on /— and d— levels. Without 
the hybridization account r and i/; i=0j ±. r are degenerate (with energy Sf < 0) and are expessed as 



?+|0>, 0+ = ^(fr+dt ~ JAY, ^±l,r = fe >r |0), ^l,r = feS"^ 



Hybridization with a singlet two / -electrons state 2j r = 2t|0), 2^ r = f£ + f£_\0) (its energy 
2£j + Uf) removes degeneration between a singlet r and triplets ipi=o t ±- r . Triplet states don't 
hybridize with singlets. As a result of hybridization between 4> r and 2f r two singlet eigen states of 
H r have the form 



Vr , ± = ^ ± \0), ft. ± = Ajt + B±2%; E v;± = ~{3e f + U f ±y/(e f + U f f + 8V 2 } 

At Uf + Efd — V we have E v = E v .- ~ Ef — V, E v - + ~ £/ + 2 V, = £/. The lower singlet 
state ip r = <^r ; - = A</>+|0) + £>2^ r |0) is a ground one. Singlet-triplet transitions will be considered at 
T < V, then the upper singlet state y? r;+ can be omitted (E^ + — E v ~ 3 V). Thus the hybridization 
role in one-site Hamiltonian if ,r is reduced to the energy renormalization of a ground singlet so that 

Etp < E^p = Ef. 

Below it appears, that instead of triplet states , 0i=o,±;r it is more convenient to use basis [IE] 

w z = tfj ; w x = - ^ = 75(^-1+^1), 

which explicitly reflects a spherical symmetry of a problem. 

In considered approximation the basis of states <p r , w qtT (q = x, y, z) is a full one-site basis. In this 
bases H (pQ) takes a diagonal form H 0s _ t an d Hj is transformed in singlet-triplet Hamiltonian H s -t 

H s -t = Ho,s-t + J ; ^o,,-t = e^Zf; (2) 

here and below Z£ x> are Hubbard projection operators into corresponding states, A, A' = ip, w q 
(Z— operators indexes q = x,y,z correspond to w q ); e — E^ — E v , below energy will counted from 
the ground state energy E v . At characteristic relation between parameters £/+ Uf — 2V the triplet 
energy is e = 0.7V, and y? r -structure (p r = A(fi r + B2f r answers the values A 2 ~ 0.8, B 2 ~ 0.2. At 
V ~ 30 meV we have for e -value e ~ 20 meV, this value is close to Ml, M2 -peaks energy. 

In the basis <p r , w q<r /—spin site operator Sj r (q = x,y,z) looks like 

Sj r = \(AZ™ + AZr + ie qq , q „Zf<i'). (3) 

In what follows in the expression (JS]) we will not distinguish the states ip r h r and will assume 
for the simplicity A = 1. 

For a two-site cluster with the Hamiltonian H s _ t ([2]) let us construct eigen states which have 
eigen values W for a cluster full spin operater. The analysis of these states allows to clarify the 
necessity of falling outside one-site approximation. The basis of four one-site operaters </>+, w^ r 
gives 16 eigen states (corresponding energies are designated as E°) with W — 0, 1, 2. 

There are two singlet states $ nim and X nm which have a form <I> n)m =0+0+ |0) > 
A n>m =^ S 9 w+ n iu+ m |0) (here |0)-cluster vucuum) with energies E% = 0, E° x = 2e. 
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Six triplet states are B\ = i(l — T n>m )u/+ n 0+ |0); D? = -^=(1 + T nim ){u+ n 0+ |0), here and below 
T n m is an operator of n, m -sites permutation. The states B\ h D\ have the same energy E B q^ D q^ = e 
but different parity T relative to T n m (T b q = — 1, Tp? = +1). 

Three more triplet states have the form B\ = ^i£qq>q"W J a wp m |0) with = 2e and T B ? = 

-1. 

At last there are five quintet states I) , I — l-r-5 with full spin W = 2, E Q Ll = 2e, parity T = (+1), 
they are described by the following wave functions 

-^(1 + T ntm )w^ n Wy m ] (1 + ^n,m)w^ n W^ m ; + Tn,m)^ n W^ m ; 

71«n^m - ^n W im)l°>; ^(^n^m + W p<m ~ 2 < n W+ m )|0>. 

The states with T = (+1) correspond to Fourier transform with k = (0), T = (—1) with k = (jr), 
the last one is an anolog of AFM point. 

The exchange J (0Q) leads to mixing between the states with coinciding T and W (and coinsiding 
q for triplets). The matrix of the Hamiltonian in $ n ,m , X nm basis gives a cluster ground state 
^ n ,m- It is formed mainly by $ n ,m which is factorable on n, m. contains an admixture of X n m 

state (with an amplitude proportional to J/e at small J), X nm state is not factorable. An explicit 
^n.m-form gives a possibility of simpliest calculations for two-site correlation functions at low T. 

The main result of the above classification is a construction of two-site nonfactorable on sites B\ 
triplet states. They are the only states that permit singlet \I/n,m triplet excitation with the energy 
E° B q = 2e. Below it is shown that exactly the movement of this two-site excitation leads to a spin 
branch which answers the position and dispersion of M3 peak. 



3 Spin susceptibility of /—electrons. 

For the description of a real spin spectrum it is necessary to turn from a two-site cluster problem 
to a consideration of an infinite lattice with the Hamiltonian if s _ 4 .((2]). The susceptibility X/(k, uj) 
is given by the two-time spin Green's function (GF) 

r 

Since YbB 12 lacks long range spin order we describe the properties of its spin system by 
a spherically symmetric approach (SSA). similar to that used in [19] . Within the limits SSA 
G^(k,u)) = Gj(k, u) and the following averages are equal to zero: (Zf q ) = 0; (Z q ' q ") = at 

Taking into account ([3]) the site representation of G z J\i, u) breaks up into a sum of terms 

(Z*\Z*) m (Z^>\Zt%), (Zt\Z?%, {Z»\Z») U , (ie zqW ,z£'«'\ie zqiq Jr 2 U (Z£ z \ie ZQiq2 Z q ^f r %, 
{ie zqlqll Z^\Zf)^ (Z^\ie zqiq2 Z™%, {ie zq , q „Z^'\Z z T *)u. 

As mentioned above the ground state is formed mainly by one-site singlets r . Then spin excita- 
tion first of all should be described by one-site singlet-triplet operator Z^r with energy e. Therefore 
below we will begin the discussion from the GF G^_ r (uj) = (Z^ Z \Z^) U which corresponds to the 
first term of full GF Gy(k, u) and has obviously the basic pole close to e. 

The equation for GF (Z^ Z \Z^) W has the form: 



[z-e)(Z** \Z*% = 5 n , r P^ + ([Zf; J\\Z?U (4) 
= P*-Pw\ P*=(Zt+); P w = {Zf). (5) 



Carrying out the commutation [Z^ z ; J] in an explicit form with the account of substitution 
of expressions ((3D in exchange interaction J, we will receive the following equation for (Z^ z \) = 
(Z^lZ^)^ (here and in obvious cases a simplified GF notation is used) 
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(z - e)(Z+*\) = S^P^ + i £ g {(i?n; g |> + [(<K;g|> - (W n;g \})+ (6) 

+2((C/ n;g |}-(K lg ;- g |})]+(5 n;g |}+(L n . g |) - (N a;s \)+2[{U+ S \) - (W+ S ._ s \) ]}. 
Here the arisen two-site operators have the form 

Un; S =Zt+ s Z^ l> ] V n . t g = Z^ s Z% d ] iy n;g = Z^Z^ d ;^i? n;g = ie Z uqZ^l^Z^ q ] (7) 

n — or 7 U< ^ 7</>Q- T — or 7 qu 7<t><i>- AT — or 7 qu 7 dz 
^n;g, ' t 2«9 Zj n+g Zj n ) n;g ' ,t 2:iig Zy n+g Zy n ) n;g l ' l ^duq Zj n+g^n " 

The GF in the right part of ((6|) describe transitions with energy uj m — e (operators , W^g), 

uj ~ (operators -D n;g , -^n ;g , ^n ;g ), w ~ e (operators f/ n;g 7 Ki;g, W'njg ) h uj ~ 2e (operator -R n; g)- 
Let us clarify their physical meaning. 

The GF with _D n;g , L n;g , iV n . g operators answer triplet-triplet against a singlet background, they 
should lead to a quasielastic peak which we donot consider. In the simplest mean field approach it 
is possible to neglect also the GF (U^. g \) , (W^-gl) as they.do not give a return to the initiale GF 

The operators t/ n;g , V n - S , W n - S (J6j) cjrrespond to the same energy transitions as Z^ z , uj = e. 

It should be noted that the term (U n . s \) describes a triplet excitation movement from site n to a 

neighbouring n + g. Operators U n - g , V n - g , W n;s permit a mean field approach with a return to GF 

(Z? z \): [((K ig |) - (Wn ig |)) +2((C/ ni g|)-(Kg;- g |))] « P^Z^Zf). In this approximation © takes 
a form 

(z - e){Z*T\) = 5 n , r P^ w + i E g {(^n;g|) + 2P+ U (Z* 9 \)} . (8) 
In a momentum representation (JSj) gives 

[z-e- jP^zMGHu) = P^ W + ^ ]T<£ k; gl>, Gi(u>) = (Zt\Z z _i) , (9) 

g 

here zi -number of the nearest neighbours, = zf 1 Eg e * kg - For the face-centered cubic YbBi 2 
lattice zi = 12; g^ = | [cos k x cos k y + cos k y cos k z + cos k z cos k x ] , the AFM vector is Q = f [1; 1; 1]. 

If to neglect the GF (i?k ; g|) in the equation ((91) it will describe a triplet excitation spectrum spin 
-^ki = e + ^P^wZigk with a gap equal e and a dispersion part ^P^ w zig^ with minimum at Q, as well 
as is observed for branch Ml. 

If in ([9]) we neglect the GF (i?k;g|) but do not restrict ourself by a mean field approximation for 
^n; g , Ki;g, ^njg -operators then the equations of motion must be introduced for the corresponding 
GF (U D - S \), (Vn ;g |), (W n;g |). We do not represent these equations because of their cumbersome form. 
Nevertheless it may be seen that their account leads to additional branches of a spectrum close to 
a pole uj ~ e. Effectively these branches correspond to Ml peak broadening with some J - scale 
structure. The last can be treated as analogue of peak M2. 

The account of a term (Ra- S \) in the equations ([6j), ([9]) is important on principle. The operator 
i? n;g = ie zuqZf^^Z^ answers to singlet-triplet transition between two-site singlet state $ nim and 
a two-site triplet state S| with energy uj = 2e. This transition lead to (Z^ Z \Z^) GF pole close to 
M3— peak energy 40meV if e ~ 20 meV. 

The equation for (-R n; g|) has a form 

(^-2e)EgA ;g l>=^ n -^ +Eg{i{2«t/ng|)_- (U n \)) - (R n J) - (M n;g |}}+ 

-W {Eb A b;g [-C'zz,g(Zn +h | ) +Cf Z ,|g-b| I )] + Eb A bjg[Cf Z, a (^n+g+b I )- C ZZ,\ S +b\ (^n+g |>]}l 

= Eg([£n; g ; ^]) ! = Ab;g = (1 ~ A b;g ) ; g = -g . 

(10) 



5 



In the equations (fTOl) right side we neglected three-site GF which give a zero input in a mean 
field approach with a return to GF , i? m;g . It is easy to see that correlation function ^ 
is a zero if it is calculated using the two- site cluster state ground state ^ n ,m- After extracting of 
averages in GF the right part (fTOl) the transition to k representation gives 

(z - 2e + |) £ g (£k; g |> = JP* Zl (l - £ k )G k (11) 
The equations ((9)) and (fTB give a final expression for 

i(z-E° ;1 )(z-E 2 )-V k 2 }G k = (z-E° 2 )Pt w , z = u + t5 
E^i = e + iP^z l9k] E\ = <le-{- V* = j^P^l - g k }. { 1 

Expression G k (fT2l) describes two triplet exitations bands i?k ; ± = HC^k-i + -^2) ^ 
\J (-^k-i + E®) 2 + 4V^}. The lower branch E k ._ (the analogue of Ml branch) is close e and has 

the dispersion with a minimum at AFM vector Q. At k = Q effective hybridization V£ (fT2]l is maxi- 
mum. Therefore the hybridization between a level E\ and E^. x gives a spectrum of the upper branch 
£k;+ (analogue of M3 branch) with a maximum at Q (as it is observed for M3 peak dispersion). 
Our consideration specify that if a system has a band in a vicinity e then with necessity it has also 
a band in a vicinity 2e. 

With temperature increase the intensity of X/(k, oj) decreases due to a factor P^ w = (Z^) — (Z qq ): 
with T -increase the expectation value (Zf^) is decreasing and (Z qq ) increases. In the limit T>e 
the expectation values (Z^), (Z qq ) tend to \ and the Ml, M3 -peaks intensity tends to zero. At the 
same time it is obviously that with the increase of the triplet state occupancy (Z qq ) the intensity of 
quasielastic peak (triplet-triplet transitions) must grow. 

Thus, our consideration allows to reflect the basic features of the dynamic magnetic response in 
YbByi [8], first of all the three-level excitation character with a specific dispersion. In the conclusion 
we will notice, that our earlier Ml peak consideration [15] is close to the present one. But it did 
not took in the account the triplet ecitation i? n;g = ie ZU gZ^1 g Z^ q , without which it is impossible to 
explain the nature of M3 peak. 

Authors are grateful to P.A.Alekseev and K.S.Nemkovski for useful discussions. 
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